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Abstract 

We  consider  the  problem  of  obtaining  integral  representation  of  feedback  operators  for  damped  hyperbolic 
control  systems.  We  show  that  for  the  wave  equation  with  Kelvin- Voigt  damping  and  non-compact  input 
operator,  the  feedback  gain  operator  is  Hilbert-Schmidt.  This  result  is  then  used  to  provide  an  explicit 
integral  representation  for  the  feedback  operator  in  terms  of  functional  gains.  Numerical  results  are  given 
to  illustrate  the  role  that  damping  plays  in  the  smoothness  of  these  gains. 


Accesion  Fey 

NTIS  CRA&!  U 

DTIC  T--.B  C 

Unannc  'need  C 

Justification 

i 

By 

Distribution/ 

Availability  Codes 

Dist 

M 

Avail  c 
Spe 

md /  or 
da! 

. 

*This  research  was  supported  in  part  by  the  Air  Force  Office  of  Scientific  Research  under  grant  F49620-93- 1-0280  and  by  the 
National  Aeronautics  and  Space  Administration  under  contract  No.  NASA-19480  while  the  author  was  a  visiting  scientist  at  the 
Institute  for  Computer  Applications  in  Science  and  Engineering,  NASA  Langley  Research  Center,  Hampton,  VA  23681-0001. 

tThis  research  was  supported  in  part  by  the  Air  Force  Office  of  Scientific  Research  under  grant  F49620-93-1-0280  while  the 
author  was  a  visiting  scientist  at  the  Air  Force  Center  for  Optimal  Design  and  Control,  Virginia  Polytechnic  Institute  and  State 
University,  Blacksburg,  VA  24061-0531,  and  by  the  National  Science  Foundation  under  grant  DMS-9409506. 


1 


1  Introduction  and  Motivation 

The  purpose  of  this  article  is  to  extend  the  representation  theorem  in  [1]  and  [7]  to  certain  classes  of  damped 
hyperbolic  systems.  The  original  motivation  for  our  study  of  hyperbolic  systems  comes  from  the  work  by 
Lupi,  Chun,  and  Turner  [8].  The  approach  in  [8]  is  interesting  because  they  make  no  prior  assumptions 
regarding  the  form  of  the  controls  and  actuators  so  that  the  gains  operators  produced  by  an  optimal  design 
could  be  used  to  make  decisions  about  where  actuators  and  sensors  are  best  placed.  In  particular,  in  [8]  it 
was  assumed  that  the  input  operator  was  the  identity  and  the  elastic  system  was  not  damped.  By  solving  the 
LQR  problem  with  the  input  operator  equal  to  the  identity,  one  can  gain  insight  into  the  type  and  location 
of  practical  distributed  controllers  for  structural  control.  This  insight  comes  from  explicit  knowledge  of  the 
kernels  (so  called  functional  gains)  that  describe  the  integral  representations  of  feedback  gain  operators. 
Even  with  no  damping  the  LQR  problem  has  a  solution  since  the  input  operator  is  the  identity  (the  system 
is  exactly  controllable).  However,  as  we  see  below  the  problem  with  no  damping  is  extremely  complex.  Basic 
questions  concerning  the  existence  and  smoothness  of  functional  gains  remain  open  and  yet  these  issues  are 
important  in  the  applications  proposed  in  [8].  Consequently,  as  a  first  step  we  take  the  middle  ground  and 
consider  damped  systems  with  distributed  control. 

Except  for  the  obvious  cases  with  bounded  input  operator  and  a  finite  number  of  controllers,  the  problem 
of  obtaining  explicit  representations  of  feedback  laws  is  more  complex  than  one  might  first  imagine.  In  the 
most  general  case,  this  problem  is  equivalent  to  the  problem  that  led  Grothendieck  to  develop  the  theory  of 
topological  tensor  products  and  nuclear  spaces.  This  theory  led  to  the  famous  Schwartz  Kernel  Theorem. 
However,  the  fact  that  the  operators  of  interest  often  arise  as  solutions  to  Riccati  equations  can  be  exploited 
to  yield  reasonable  results.  The  representation  problem  for  a  parabolic  problem  with  unbounded  control 
operator  was  first  considered  in  [1]  and  [7].  It  was  shown  that  as  long  as  the  input  operator  is  bounded 
relative  to  the  open  loop  dynamic  operator,  the  solution  to  the  algebraic  Riccati  equation  is  Hilbert-Schmidt. 
This  fact  was  then  exploited  to  show  that  the  resulting  feedback  operator  had  an  integral  representation. 
For  the  hyperbolic  case  considered  here,  these  types  of  results  are  more  intricate  and  highly  dependent  on 
the  type  of  damping. 

In  the  Section  2  we  present  a  numerical  example  involving  the  control  of  a  hybrid  cable-mass  system. 
This  example  is  used  to  motivate  the  model  problem  and  to  demonstrate  that  distributed  controllers  can 
enhance  disturbance  attenuation.  We  then  concentrate  on  the  ID  wave  equation  with  Kelvin-Voigt  damping. 
Although  similar  results  can  be  obtained  for  other  damping  models,  we  present  the  basic  theorem  for  this 
model  in  order  to  keep  this  paper  moderately  short.  However,  we  illustrate  the  ideas  and  difficulties  for 
other  damping  models  with  numerical  examples.  Finally  we  close  with  a  few  comments  about  future  work 
and  other  open  problems. 

2  Control  for  a  Cable-Mass  System 

The  following  system  was  proposed  by  Nayfeh,  Nayfeh  and  Mook  [11]  as  a  simple  example  of  a  nonlinear 
distributed  parameter  system  with  the  property  that  many  standard  discretized  lumped  models  failed  to 
capture  the  essential  nonlinear  behavior  of  the  dynamic  system  governed  by  the  partial  differential  equation. 
This  system  was  also  considered  in  [2]  where  it  was  used  as  a  test  model  for  M INMAX  control  for  distur¬ 
bance  attenuation.  However,  in  [2]  there  was  only  one  controller  and,  although  disturbance  attenuation 
was  achieved,  we  shall  see  below  adding  a  distributed  controller  can  improve  performance.  The  MINMAX 
approach  provides  a  “robust  state  feedback  control  law”  which  is  less  sensitive  to  disturbances  and  certain 
unmodeled  dynamics  than  is  the  LQR  design.  The  idea  is  to  obtain  a  representation  of  the  control  law 
and  then  use  approximation  theory  to  compute  finite  dimensional  sub  optimal  controllers.  These  suboptimal 
controllers  were  used  to  attack  the  problem  of  designing  reduced  order  state  estimators.  We  shall  limit  our 
discussion  here  to  the  full  state  feedback  problem. 

Consider  the  hybrid  nonlinear  distributed  parameter  system  described  by  a  vibrating  cable  held  fixed  at 
one  end  and  with  a  mass  attached  at  the  other  end.  The  mass  is  suspended  by  a  spring  which  has  nonlinear 
stiffening  terms  and  is  forced  by  a  disturbance  (see  Figure  2.1).  The  equations  for  the  hybrid  system  are  [2] 


1 


(2.1) 

(2.2) 


The  initial  conditions  are  given  by 


w(t,  0)  =  0. 


(2.3) 


w(0,s)  =  w0,  —w(0,s)  =  w\ . 


Here,  w(t,s)  represents  the  displacement  of  the  cable  at  time  t,  and  position  s ,  w(t,£)  represents  the 
position  of  the  mass  at  time  t,  p  and  m  are  the  densities  of  the  cable  and  mass  respectively,  r  is  the  tension 
in  the  cable,  and  7  is  a  damping  coefficient.  The  alphas  are  coefficients  describing  the  nonlinear  effects  of  the 
spring.  The  term  rj(t)  is  viewed  as  a  disturbance  and  u\(t,  s )  and  U2(t)  are  control  inputs.  For  the  moment 
we  assume  Kelvin- Voigt  (internal)  damping  in  the  cable. 

The  problem  is  hybrid  in  that  the  system  is  described  by  a  linear  partial  differential  equation  (the  wave 
equation)  coupled  through  the  boundary  condition  to  a  low  order  nonlinear  ordinary  differential  equation 
(Duffing’s  equation).  In  [2]  it  was  assumed  that  the  control  acted  exclusively  on  the  mass  (i.e.  that  ui(t,x) 

=°). 

This  model  is  often  first  written  as  a  second  order  system  in  a  Hilbert  space  H  of  the  form 


Figure  2.1:  Cable-Mass  System 

This  formal  system  has  the  advantage  that  it  has  the  same  appearance  as  the  finite  dimensional  case 
and  in  order  to  address  viscous  and  “structural”  damping  one  merely  replaces  Do  =  Ao  with  D0  =  I  and 
Do  —  [A0]lf2,  respectively.  We  note  however,  that  it  is  more  consistent  with  physics  to  write  the  system  in 
the  form 

y(t)  +  S*(Sy(t)  +  7  Ty(t))  +  *b(y(0)  =  *o«(i)  (2.6) 

where  S  =  T  =  [^o]1^2*  Observe  that  Aq  =  Aq  >  0  and  so  S*  =  S  and  S*S  =  S*T  =  Ao>  Hence,  (2.6)  is 
formally  obtained  by  factoring  [^o]1^2  out  of  the  expression  D0y(t)  +  A0y(t)  in  (2.5).  Note  also  that  (2.6)  is 
of  a  form  that  allows  for  structural  damping  where  S  =  [Aq]1/2  and  T  =  I,  as  well  as  for  viscous  damping 
where  5  =  [^o]1^2  and  T  =  5_1  =  [ Ao J”1/2.  In  addition,  by  writing  the  system  in  the  second  order  form 
(2.6),  one  captures  a  form  that  comes  from  balance  laws  and  at  the  same  time  sets  the  stage  for  a  simple 
formulation  of  the  problem  in  first  order  state  space  form. 
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The  system  governed  by  equations  (2.1)  -  (2.4)  can  be  written  as  a  dynamical  system  in  an  appropriate 
(infinite  dimensional)  state  space.  Although  there  are  several  equivalent  formulations  for  this  problem,  we 
shall  write  the  governing  equations  as  the  first  order  system 

x(t)  =  Ax(t)  +  F(x(t))  +  Bu(t)  H»  Dr](t),  x(0)  =  x0  (2.7) 

where  at  time  t  the  state  x(t)  =  [y(t),  y(t)]T  lies  in  the  Hilbert  space  X  =  H l  x  1R  x  L2  x  IR.  Here,  Hi  is  the 
subspace  of  the  Sobolev  space  H1  =  ^ 1  [0 ,  defined  by  H[  =  {w  £  H1  :  w( 0)  =  0},  and  L 2  is  the  standard 
Lebesgue  space  of  square  integrable  functions.  The  control  u(t )  lies  in  the  control  space  U  =  L2  x  IR.  Here 
the  inner  product  in  X  is 


[*(•). I =  T  f  w'(x)w'(x)dx  +  aitt 

JO 

Ap  /  vt{x)v\x)dx 

Jo 


ix  -f*  mptjji. 


(2.8) 


It  is  important  to  precisely  define  the  system  operators  and  their  domains  in  order  to  obtain  correct  rep¬ 
resentations  of  the  feedback  operators  that  will  be  used  to  control  the  system.  Let  61  denote  the  “evaluation 
operator”  defined  on  H 1[0,i]  by  ^(<£(*))  =  </>(£)  and  define  the  linear  operator  A  on  the  domain  V{A)  C  X 
by 


2?(J4)=( x  =  [w,£,v,n]T  eX:w,ve  Hi,  1“^^+ 

w(t)=Z,v(£)=  (i), 


and 


[  d  f  r  d  7  d  }  c  f  r  d  )  d  ) 

A’=  [v^ts  \fTsw+ ?**/  -*1  ;r+  m*”r 


m 


The  control  input  operator  B  and  the  disturbance  operator  D  are  defined  by 


B  =  [0,Ih]  and  Dr]  = 
respectively.  The  nonlinear  operator  F  is  defined  on  X  by 


°,M,—  r] 

m 


F(x)=\0,0,Q-^[e]T  =  [0,Fo(y)f. 


(2.9) 


(2.10) 


(2.11) 


(2.12) 


Observe  that  the  input  operator  is  the  same  as  that  used  in  [8]  in  their  analysis. 

As  noted  above,  this  problem  with  ui(<,$)  =  0  was  considered  in  [2]  where  MINMAX  control  was  used 
to  design  a  low  order  dynamics  control  law.  This  law  was  based  on  two  outputs  (position  and  velocity  of 
the  mass)  and  resulted  in  a  practical  low  order  design.  For  this  note,  we  shall  consider  only  the  full  state 
feedback  problem.  However,  we  allow  for  distributed  control  through  ui(£,s).  The  simplest  approach  is  to 
linearize  the  system,  use  MINMAX  design  to  obtain  a  feedback  operator  and  apply  this  law  to  the  full  plant. 

The  linearized  system  has  the  form 


— x(t)  =  Ax(t)  +  Bu(t)  +  Dr)(t ),  x(0)  =  xQ. 


(2.13) 


For  this  problem,  with  Kelvin- Voigt  damping,  one  can  apply  the  MINMAX  theory  in  [10]  to  obtain  a  feedback 
law  of  the  form 


= 


Ui(t,s) 
«2  (0 


-Kex(t) 


(2.14) 


where  for  9  >  0  the  gain  operator  Kg  :  X  — *  H  is  given  by 


Ke  =  B*Pe 


(2.15) 
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(2.16) 


and  Pq  satisfies  the  algebraic  Riccati  equation 

A*P  +  PA-  P[BB*  -  02M]P  +  Q  =  0. 

Here  M  =  M*  >  0,  Q  =  Q*  >  0  and  the  system  (2.16)  holds  in  the  weak  sense  (see  [10]).  When  0  =  0 
one  has  the  LQR  design.  However,  when  0  >  0  the  corresponding  MINMAX  controller  provides  additional 
disturbance  attenuation  (see  [2]). 

Figure  2.2  clearly  illustrates  the  difference  between  the  performance  achieved  in  [2]  and  what  one  can 
achieve  with  additional  distributed  control.  For  this  example,  we  use  the  same  parameters  and  finite  element 
scheme  found  in  [2] . 
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When  there  is  control  on  the  mass  only,  the  maximum  value  of  0  yielding  a  feedback  law  of  the  form  (2.15) 
-  (2.16)  was  0  =  1.7;  this  behavior  is  shown  by  the  dashed  line.  If  in  addition,  one  allows  distributed  control 
(at  0  =  1.7)  there  is  a  loss  in  attenuation  (but  possibly  better  performance)  as  shown  by  the  dotted  line. 
However,  by  allowing  control  on  the  cable,  the  value  of  theta  can  he  increased  to  0  =  2.5,  greatly  increasing 
disturbance  attenuation  and  improving  performance  as  shown  by  the  solid  line  in  Figure  2.2. 


Figure  2.2:  Displacement  of  the  Mass  under  Three  Types  of  Control. 


This  example  clearly  indicates  the  benefit  of  distributed  control  and  MINMAX  design.  Moreover,  as 
shown  in  [2]  if  one  has  explicit  representation  of  the  feedback  operator  K$,  then  this  can  be  used  to  design 
practical  low  order  dynamic  controllers  (nonlinear  compensators).  In  particular,  it  follows  that  since 


'  K]1  Kg2  Kf  Kj 4 
Kf  Kf  Kf  Kj 4 


is  bounded  from  X  to  H,  the  control  on  the  mass  has  the  form 

3 

u2 (t)  =  -  J  rks(z)-^w(t,z)dz  -  aika(t,£) 

3  3 

pkv(z)-j^w(t,z)dz  -  mkm-^w(t,£). 


(2.17) 


(2.18) 
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It  is  tempting  to  assume  that  one  also  has  a  representation  of  the  form 


fl  d 

(t,C)  =  ~  J  rhs(C,  z)dz  -  aiha(()(t, 

d  d 

phv(C,  z)—w(t,  z)dz  -  mhm(C)—w(t,£). 

In  (2.18)  k$(z)  and  kv(z)  are  functional  gains  corresponding  to  strain  and  velocity,  respectively.  Their 
existence  and  smoothness  properties  are  assured  by  the  Riesz  Representation  Theorem.  For  example,  there 
is  a  ks(z)  e  L2( 0,  £)  so  that  iff1  :  Hi  — ►  IR  has  the  form 

[■ K I1]  <£(■)  =  rks(z)-j-<f>(z)dz.  (2.20) 

However,  on  the  surface  all  we  know  about  iff3  :  L2(0, 1)  — ►  1/2(0, 1)  is  that  it  is  bounded.  As  noted 
above,  the  desire  to  find  representations  of  such  operators  as  integrals  led  to  Grothendieck’s  work  on  nuclear 
spaces  and  the  Schwartz  Kernel  Theorem.  Recall  that  not  all  bounded  linear  operators  on  1/2(0, 1)  (even  if 
self-adjoint)  have  integral  representations,  as  illustrated  by  the  identity  operator. 

In  [8]  this  issue  was  avoided  by  assuming  a  representation  similar  to  (2.19)  and  then  allowing  generalized 
functions  as  kernels.  This  approach  proved  to  be  satisfactory  for  the  one  dimensional  case  considered  therein, 
but  does  not  apply  to  more  general  two  and  three  dimensional  hyperbolic  problems.  We  present  some  results 
that  lead  to  a  representation  of  the  form  (2.19)  when  there  is  suitable  damping  in  the  system.  Although  this 
approach  does  not  apply  to  the  undamped  case,  the  ideas  can  be  extended  to  certain  higher  dimensional 
damped  elastic  systems.  Moreover,  we  conjecture  that  the  undamped  problem  in  3D  systems  will  not  yield 
a  representation  even  of  the  type  considered  in  [8]. 

We  turn  now  to  the  simple  ID  wave  equation  in  order  to  state  precise  results  and  to  keep  this  article 
at  a  reasonable  length.  Moreover,  we  restrict  our  presentation  to  LQR  design  (0  =  0)  since  we  can  rely  on 
existing  literature  to  outline  the  results.  Extensions  to  higher  dimensional  problems  with  0  >  0  will  appear 
in  a  future  paper. 


(2.19) 


3  The  Wave  Equation 


Consider  the  wave  equation  with  Kelvin- Voigt  damping.  Damping  plays  a  key  role  in  the  design  of  controllers 
for  hyperbolic  systems.  To  illustrate  this  role,  consider  the  LQR  control  problem  for  the  wave  equation 
defined  by 

6 2  d 2  d 3 

$ )  =  S)  +  idsidi^’ s)  +  u(<) 


with  boundary  conditions 
and  cost  function 


0  <  s  <  1,  0  <  t, 
w(t,  0)  =  0,  w(t,  1)  =  0,  0<t 


1  r°°  r1  d 

■'w  =  u  1 


+ 


9  u  \ 
— uj(<,  s) 


+  Mm)P 


dsdt. 


(3.21) 

(3.22) 

(3.23) 


This  problem  is  defined  on  the  state  space  X  —  #o(0, 1)  x  L2(0 , 1).  As  above,  care  must  be  used  to  define 
the  system  operators.  In  particular,  let 


V(A)  =  (  wv  €A-:»>w€^(0,l)>ti;  +  ‘ri>€ir1(0>l)| 


and  define  A  by 


w 

v(-) 

V 

.  Ts  [-hw  +  'lTsv\  . 

(3.24) 

(3.25) 
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The  control  space  is  U  =  £2(0, 1)  an<i  the  control  operator  is  defined  from  U  into  X  =  #0  x  £2  by 

5  =  [  r°  1  .  (3.26) 

■^£'2 

Here,  0  =  0  and  Q  =  R=  Ix.  The  LQR  problem  has  a  solution  (even  if  7  =  0)  given  by 

■»•■)—*[*&?.)]•  <3'27) 

where  K  :  i/g(0, 1)  x  1)  has  the  form 

K  =  B*P  (3.28) 

and  P  satisfies  the  weak  form  of  the  Riccati  equation  (ARE)  given  by 

(Par,  Az)x  +  (Aar,  Pz)x  -  (P*Pz,  P*P*)c/  +  (Car,  C*)y  =  0,  (3.29) 

for  all  x,y  in  2>(A).  For  this  second  order  system,  K  takes  the  form 

K  =  [0,  ILa]P,  (3*30) 


or  equivalently, 

K  =  [P2 u  P 22]  (3.31) 

where 

P  = 

The  operators  P2i  :  if^(0, 1)  £2(0, 1)  and  P22  :  T2(0, 1)  ->  L2( 0, 1)  are  bounded  linear  operators  with 

P*2  =  P22.  The  goal  is  to  determine  if  there  exist  “nice”  integral  representations  of  these  operators. 

The  following  result  is  well  known  and  may  be  found  in  [3] . 

Lemma  3.1  The  operator  A  generates  an  analytic  semigroup  S(t)  on  X  and  there  exist  M  >  0,u>  >  0  with 
||S(i)||  <  Me~wt. 

Since  A  generates  a  stable  analytic  semigoup,  we  can  apply  Theorem  2.1  (page  36)  in  [9]  to  obtain  the 
following  regularity  result  for  P. 

Lemma  3.2  There  exists  a  self-adjoint7  non-negative  definite  bounded  linear  transformation  P  =  P*  satis¬ 
fying  (3.29).  Moreover ,  for  each  e  >  0,  the  operators  [A*]1“€P  belong  to  C(X,X). 

We  note  that  Theorem  2.1  in  [9]  also  states  that  e  can  be  set  equal  to  zero  if  A  is  self-adjoint,  normal,  or 
has  a  Riesz  basis  of  eigenvectors.  However,  it  is  interesting  to  note  that  A  is  neither  normal  nor  self-adjoint. 
For  this  ID  problem,  A  does  have  a  Riesz  basis.  However,  this  property  is  not  needed  to  establish  the 
following  representation. 

Theorem  3.3  There  exist  functions  &2i('>  *)>  ^22(*>  *)  such  that 

(1)  M-,  •)  €  L2([ 0, 1]  x  [0, 1]),  M->  •)  €  L2([ 0, 1]  x  [0, 1]), 

(2)  k22(C,s)  =  k22(s,C), 

(3)  the  mapping  t  — ►  k2 ,(£,  t)  belongs  to  Hi  for  almost  all  £  6  [0, 1] 
and  one  has  the  representations 

k21(C,t)<f>(t)dt,  <f>eHl,  (3.33) 

[P22<t>]( 0=  f1  k22((,t)<f>(t)dt,  <t>  €  L2.  (3.34) 

Jo 


[P21<f>](  C)=  [ 

Jo 


Pll  Pl2 
P21  P22 
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The  proof  of  Theorem  3.3  is  rather  tedious  and  will  not  be  included  here.  However,  we  note  that  the 
proof  is  similar  to  the  proof  of  Theorem  3.2  given  in  [7].  It  relies  on  a  classical  theorem  by  Fullerton  (see 
Theorem  6  in  [5]).  However,  the  specific  structure  of  the  A  operator  combined  with  the  special  B  operator, 
B  =  [0  Ijr]T,  is  needed  to  carry  out  the  proof  for  this  hyperbolic  case.  Theorem  3.3  leads  easily  to  the 
following  representation. 

Theorem  3.4  There  exist  functional  gains  kv( *,  •)  and  £*(-,  *) 

(1)  0  €  c([ o,  i]  X  [0,  i]),  *,(-.  •)  e  MO,  i]  x  [o,  i]), 

(2)  kv(C,z)  =  kv(z,  C) 

(3)  the  mapping  z  — ►  &$(£,£)  belongs  to  H2  for  almost  all  £  E  [0, 1] 
and  the  LQR  feedback  control  law  has  the  representation 


fl  d  f1  d 

u(t,C )  =  -y  ks(C,z)-w(t,z)dz-J ^  kv(C,z)-^w{t,z)dz 


Proof:  Since  K  :  Hq  x  L%  —*  L?  is  given  by 


’(■)  =~Jo  k21^’z^w^dz~  J0  k^’z)v(z)ds 


let  •)  =  Ar22(*,  ■)■  Integration  by  parts  on  the  first  integral  yields  (for  w( ♦)  E  Hq) 


[  Ar2i(C ,z)w(z)dz=-  [  \[  k2i((yt)dt  w\z)dz . 
Jo  Jo  LJo 


The  representation  (3.35)  follows,  where 


ks(C>z)  =  -  [  hl(Ct)dt, 
Jo 


(3.35) 


(3.36) 


and  2  ks(C,  z)  belongs  to  H2. 

We  turn  now  to  some  numerical  experiments  to  illustrate  the  representation  (3.35).  In  addition,  we 
consider  other  damping  models  and  present  numerical  results  that  clearly  show  the  role  that  damping  plays 
in  the  smoothness  and  existence  of  functional  gains. 


4  Numerical  Experiments 

We  consider  the  wave  equation  with  various  forms  of  damping.  We  have  the  representation  (3.35)  for  Kelvin- 
Voigt  damping  and  similar  results  can  be  obtained  for  structural  damping.  However,  as  we  shall  see  below, 
viscous  damping  is  not  sufficient  to  ensure  the  existence  of  L2  functional  gains  even  though  the  feedback 
operator  exists  and  is  bounded.  We  set  7  =  .25  and  use  standard  piecewise  linear  finite  elements  to  compute 
ks( •,*)  and  kv (-,*).  In  particular,  we  use  finite  elements  to  approximate  the  Riccati  equation  (3.29)  and 
construct  K  as  defined  by  (3.30-3.32). 

In  order  to  show  convergence  of  the  scheme  we  compute  A;^(-,  •)  and  £^(-,  *)  for  Kelvin- Voigt  damping. 
Here,  N  represents  the  number  of  elements  in  the  model  (i.e.,  [0,1]  is  partitioned  into  (N-bl)  subintervals). 
Figure  4.1  shows  the  convergence  of  fc^(*>  *)  an<^  ^v(mi  0  to  M*>  ’)  and  M*>  *)>  respectively.  Observe  that  the 
convergence  of  fe^(*,  •)  is  very  rapid. 

Although  we  presented  theoretical  results  for  the  case  of  Kelvin- Voigt  damping  only,  Figure  4.2  shows 
that  similar  results  hold  for  structural  damping.  In  particular,  an  integral  representation  exists  and  the 
finite  element  approximations  converge.  It  is  important  to  note  that  the  functional  gain  for  strain,  k$( *,  •), 
remains  smooth  and  is  the  same  as  the  gain  obtained  with  Kelvin- Voigt  damping.  However,  there  is  a 
marked  difference  in  the  smoothness  of  fcv(-,  •).  This  functional  gain  has  a  sharper  “peak”  at  £  =  z  than 
the  corresponding  velocity  gain  for  Kelvin- Voigt  damping.  Although  we  expect  that  for  structural  damping, 
&v(C>  *)  belongs  to  H1,  we  conjecture  that  t  — »  kv(Cyt)  is  not  H2 . 


Finally,  for  viscous  damping,  the  numerical  results  shown  in  Figure  4.3  show  that  again,  k ^ (*,  •)  is  well 
behaved  and  yet  k^(-,  •)  seems  to  be  as  singular  measure  concentrated  at  (  =  z.  Thus,  we  conjecture  that 
the  representation  (3.35)  does  not  hold  for  the  wave  equation  with  viscous  damping  for  any  L2  function 
kv( •,*).  Similar  results  were  noted  in  [8]  for  undamped  beam  equations  and  in  [7]  for  parabolic  equations 
with  highly  unbounded  input  operators. 

The  numerical  results  also  indicate  that  the  strain  functional  gains  are  independent  of  the  damping  model. 
As  shown  in  Figure  4.4,  k^(*,  ♦)  does  not  change  as  the  damping  model  changes.  Thus,  it  seems  from  the 
numerical  experiments  that  damping  has  the  most  impact  on  the  existence  and  smoothness  of  kv{ •). 

5  Conclusions 

In  this  paper  we  have  provided  an  integral  representation  theorem  for  the  LQR  feedback  operator  for  the  ID 
wave  equation  with  Kelvin- Voigt  damping  and  fully  distributed  control.  The  motivation  for  this  effort  comes 
from  the  problem  of  designing  suboptimal  low  order  dynamic  compensators  and  for  optimal  sensor  location. 
Although  the  theoretical  results  presented  here  are  restricted  to  problems  with  Kelvin-Voigt  damping,  the 
numerical  results  suggest  that  more  general  results  are  available.  Also,  these  same  numerical  results  indicate 
that  once  the  damping  is  insufficient  to  ensure  the  analyticity  of  the  control  system,  the  existence  of  L2 
functional  gains  is  in  doubt.  In  particular,  the  solution  of  the  Riccati  equation  may  not  be  Hilbert-Schmidt. 

The  LQR  problem  for  distributed  parameter  systems  has  received  considerable  attention  during  the  past 
ten  years.  However,  problems  in  which  B,  Q  and  R  are  all  non-compact  have  not  been  fully  explored.  Temam 
[12]  considers  the  differential  Riccati  equation  with  B  =  R  =  7,  but  assumes  that  Q  is  Hilbert-Schmidt. 
Similar  results  are  found  in  [4,  6].  Finally,  as  the  numerical  evidence  suggests,  many  theoretical  issues  are 
not  yet  settled. 
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k_s(.,.),  N  =  8  k_v(.,.),N  =  8 


Figure  5.1:  Convergence  of  the  Functional  Gains,  Kelvin- Voigt  Damping. 
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k_s(.,.),  N=8,  Kelvin-Voigt 


k_s(.,.),  N=16,  Kelvin-Voigt 


0  0 


k_s(.,.),  N=8,  Viscous 


k_s(.,.),  N=16,  Viscous 


Figure  5.4:  Structural  Gains  for  Various  Damping  Types,  N  =  8,16 
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